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Percolation on finitely ramified lattices, those where the loss of a finite number of bonds can 
forestall end-to-end transport when the number of nodes N increases, is generally considered to be 
trivial^. They percolate only at full bond density, p = 1. Examples are the one-dimensional lattice 
and fractals such as the Sierpinski gaskeli^. We show here that by dressing up such lattices with 
small-world bonds-, in a way that the resulting network becomes infinitely ramified, an entirely 
novel percolation transition at a critical point pc < 1 emerges. Approaching this transition from 
below, many clusters of already diverging, yet, sub-extensive, sizes combine explosively^ at pc such 
that the resulting cluster becomes extensive. That is, the usual order parameter Poo, describing 
the probability of any node to be part of the largest cluster, jumps instantly to a finite value at 
Pc- Simple examples of this transition are provided by small- world networks consisting of a one- 
dimensional lattice combined with various hierarchies of long-range bond*"! that reveal many features 
of the transition explicitly. 



The percolation properties of networks are of sig- 
nificant interest; without percolation, any large-scale 
transport or communication through the network ceases. 
Much research has been dedicated to the understanding 
of percolation on randomly grown, complex networks^. 
Yet, the engineering of artificial networks with well- 
controlled features seems desirable. Indeed, there has 
been considerable interest in the properties of spatial 
networks, linking real-world geometry with small-world 
effects^i^. In particular, networks possessing hyperbolic 
and/or hierarchical features^iiiSrii relate to actual trans- 
port systems such as for air travel, routers, and social 
interactions. Hierarchical networks exhibit a number of 
novel features, as shown with exact methods below, most 
notably a discontinuous transition in the formation of an 
extensive percolating cluster for ordinary, random bond 
addition. For random networks, such an "explosive" per- 
colation transition, one that requires correlated bond ad- 
ditions, has been proposed by Achlioptas et al.^ only re- 
cently. Intense research is ongoing into the dynamics of 
correlated cluster formationi^""— , with the existence of 
a discontinuity still in doubt^^"— . Understanding of ex- 
plosive cluster formation could be significantly advanced 
when the discontinuity can be studied rigorously, merely 
by adding bonds randomly on an engineered network. 

A simple example of a hierarchical network is devel- 
oped in Figllja). In its recursive construction, bonds 
from the preceding generation are built up in the next 
generation according to a fixed pattern up to a gener- 
ation fc — >■ oo. Applying a venerable technique called 
real-space renormalization group (RG)'^^ to the complete 
network in effect reverses the recursive built-up of the 
network, as shown in Figdlb). Starting with a network 
of size N = 2'^ -1-1, the probability of two neighboring sites 
to be directly linked is Tq — p. From that, we can recur- 
sively determine the probability T„, < n < fc, that two 
sites originally separated by 2" positions in the horizon- 



tal direction have a direct path between them. Finally, 
Tn=k provides the probability for any end-to-end path in 
the network; for the two end-points all paths are direct. 
As explained in Figl^b), these probabilities satisfy the 



recursion 



p+{l-p)Tl in^p). 



(1) 



For infinitely large networks we find limfe_^oo Tn=k = T* 
from the stationary points, e.g., Eq. ([T|) yields the fixed 
point T* = 1 for all p and a line of fixed points. 
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that intersect at p = i. Analysis shows^ that Eq. ^ 
is a stable fixed point for all < p < ^ which describes 
the asymptotic behavior of large networks with a finite 
probability < T* (p) < 1 to possess an end-to-end path. 
For i < p < 1, there is always such a path, i.e., T* — 1. 

It is remarkable that percolation on these networks 
has a number of features entirely absent from basic 
lattices^ii^ii^. Specifically, on a lattice (like the Id line) 
the existence of an end-to-end path is synonymous with 
the existence of an extensive cluster, giving each node a 
finite probability Poo (p) to be connected to it. This is 
not generally true for these hyperbolic-shaped networks. 
Typically, there are two transitions: below the lower 
transition pi there are only finite clusters and no end-to- 
end paths, and above the upper transition j)„ > pi there 
is an extensive cluster and a definite end-to-end path, 
as on any lattice. However, between pi and p„ exists 
a regime with diverging but sub-extensive clusters and 
a variable probability for an end-to-end path, as given 
by Eq. Those sub-extensive clusters scale with a 

new fractal exponent < 5* (p) < 1.'^*' Accordingly, the 
order parameter Po^iil is strictly non-zero only above 
Pit, making it the thermodynamically correct transition 
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Figure 1: (a) Pattern for the recursive generation of a hier- 
archical network (displayed at full bond density, p = 1). In 
k = 0-th generation, the network consists of a single bond 
between two end-nodes (open circles). In successive genera- 
tions k — 1,2, . . ., each (horizontal) lattice bond is replaced 
by two of its kind, with a new (black) node in-between, and a 
long-range bond (shaded arcs). In generation k the network 
has A*' = 2* -I- 1 nodes and L = 2'°'^^ — 1 bonds for a (sparse) 
average degree 2L/N ^ 4 for large k. (b) The n-th RG step 
reduces each three-node graph-let into a single lattice bond 
of renormalized percolation probability for a direct path be- 
tween nodes r„+i, which derives from that probability, r„, 
of the prior bonds and the probability p to have a long-range 
bond. r„+i percolates if either the long-range bond exists 
(probability p, irrespective of Tk) or both lattice bonds are 
present (probability (1 -—p)T^), leading to Eq. ([T]). Revers- 
ing the pattern set in (a), the RG starts at n = from an 
unrenormalized network of generation k with uniform prob- 
abilities p for the long-range bonds as well as for the direct 
path to neighboring nodes, To = p. 




Figure 2: (a) Schematic for generating functions Tn{x) and 
Sn{x,y). The open circles represent the end-nodes, shaded 
areas indicate clusters that either span (T„) or do not span 
(Sn) between the end-nodes. Clusters that do not reach 
an end-node are ignored, (b) All diagrams contributing to 
T„+i{x) or Sn+i{x,y) in the nth RG step. The remaining 
end-nodes (always-open circles) are not counted in the gener- 
ating functions; the (black) connecting nodes increment the 
cluster size and are taken into account by multiplying x or 
y. The contribution of each graph-let is: (i) xpT^{x), (ii) 
xpT„(x)S„{x,x), (in) xpT„{x)Sn{x,x), (iv) pS„{x,l)S„(l,x), 
(v) x(l-p)T^{x), (vi) x{l - p)T„{x)S„{x,y), (vii) y{l - 
p)T„{y)S„{x,y), and (viii) (1 - p)S„(a;, 1)S„(1, y). As tal- 
lied up in Eqs. (|5I6|) . graph- lets (i)-(v) span end-to-end and 
contribute to Tn+i{x), while (vi)-(viii) do not span and con- 
tribute to Sri+i{x,y). 

point, pc = Pu- In the present network it is pi ~ 0; a 
more elaborate networl i^'^^ with non-trivial pi and p^ is 
mentioned below. 

To reveal the nature of the transition, we obtain the 
exact Poo from cluster generating functions. To this end, 
we introduce two basic quantities: the probability t^f'^p) 
that both end-nodes are connected to the same cluster 
of size i, and the probability s|"-'(p) that the left (right) 
end- node is connected to a cluster of size i (j) but these 
clusters are not the same. For convenience the end-nodes 
are not counted in the cluster size, as shown in Figl^Ja). 
The corresponding generating functions are defined as 

00 

T„(x) = (3) 

00 00 

Sn{x,y) = Y.J2tj(pyy'- (4) 

The recursion relations for these generating functions 
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Figure 3: Plot of Poo (p), evaluated after n = 10^ — 10* iter- 
ations of the recursions in Eq. (|10|) . evolving into a discon- 
tinuity at p — J> pc = I with Poo (Pc) = 0.609793 .... Con- 
vergence is slowest just below pc, since finite-size corrections 
decay slowly with Ar*(p)-i with 1 - * (p) ~ {pc - p) / In 2 for 
p ^ Pc from Eq. ([T3)) . 



can be obtained by considering all possible triangular 
graph-lets, like in Fig. [TJb), but now also taking clus- 
ter sizes into account, as shown in Figl^Jb). For each 
graph- let, the type of connectivity between each of the 
nodes is checked, see Figl^^a), and the bonds are assigned 
a value Tn{x) or Sn{x, y), depending on the type of clus- 
ter each represents. As in Fig. [TJb), small world bonds 
are merely assigned the probability p or 1 — p for being 
present or not. Marking the increment in cluster size, 
the inner node provides a factor of a; or y for its adjacent 
cluster, or unity if it remains isolated. The contribution 
of each graph-let to the next generation is the product 
of the weights of the three bonds and the intermediate 
node; all eight of these are determined in FiglJJb). These 
provide the recursions 

Tn+i{x) = xT^{x) +p[2xTn{x)Sn{x,x) 

+5„(x,l)5„(l,x)], (5) 
Sn+iix,y) = {1 - p)[xTn{x)Sn{x,y) + 

+yTn{y)Sn{x,y) + Sn{x,l)Sn{l,y)] (6) 

initiated with Tq(x) — p, So{x,y) — 1 — p. Note that 
both equations reduce to Eq. ([!]) for x = y — 1, where 

Eqs. dSH) give r„ (1) = 1 - Sn (1, 1) - r„. 

It is convenient to rewrite the recursions in Eqs. (|5l6p 
in terms of functions of a single variable x only. To 
this end, we define the functions S„(a;) = S'„ (x^x) and 
Sn{x) = Sn (x, 1) — Sn (1,2;) that, combined into a vec- 
tor F = (T, lead to 



Vn 



_i (x) = f (k (x) , x) 
with a function-vector T of non-linear components 



J" 



xT^ + 2xpTJ: + pS^ 
(v,x) = (1-p)(2.tTS + 52) 
^ ^ * il-p){l + xT)S 



(7) 



(8) 



The mean size (Smax)^. of the (likely largest) cluster 
connected to the end-nodes and, hence, the order param- 
eter Poo{p) = hmfc^oo (smax)fc/A^fc result from the first 
momeirt of the geirerating functions via the first deriva- 
tive in x at a: = 1, 



= T!+ Si, ^ TV, 



*(p) 



(9) 



for a network of size A^^ = 2*^ + 1 ^ oo. Expanding Eq. 
([7]) to first order in e = 1 — x — > yields to zeroth order 
Eq. ^ evaluated at x = 1. such that each component 
simply reproduces Eq. ([T]) again. To order e, we find an 
inhomogeneous linear recursion for = {x = 1) 



v;;+, 

with the Jacobian matrix 



dV 



2T + 2pY., 2pT, 2pS 

2(1 -p)S, 2(1 -p)r, 2[l -p)S 

(l-p)S, 0, (l-p)(l + T) 

and inhomogeneity 



5a; V' 



T'^ + 2pTS \ 
2(l-p)ri] 
{l-p)TS ) 



(12) 



which depend on n through V^„ (a; = 1) = 
(r„, 1 — T!„, 1 — T„). Since each network at n 
only consists of end-nodes, all clusters are empty and 

n = 0. 

For large n, i.e., near the fixed point V* = 
(T*, 1 - T*, 1 - T*), it is easy to show that the inho- 
mogeneity in Eq. (jlOp is subdominant, leaving a simple, 
homogeneous, linear system. For each component of V|^^ 
the dominant contribution then arises from the eigen- 

state with the largest eigeirvalue of ll (v*,l). Below 

the transition, < p < Pc = |, we use Eq. ([2)) for T* 
to obtain as the largest eigenvalue A = 1 -I- 2p, such that 
(smax)fc ^ A''" for fc — > oo. Then Eq. ^ yields the fractal 
exponent 



(13) 



for sub-extensive cluster sizes, with < ^ < 1, implying 
Poo = here. 

Above and at the transition, Pc < P ^ 1, = 1 and 
Eq. (fTTj) provides uniformly A = 2 as the largest eigen- 
value (i.e., = 1), indicating the existence of an exten- 
sive cluster. But unlike for a continuous transition, where 
-Poo [p) — ^ for p — 7> p+, Eq. pn]) provides for a mono- 
tone increasing probability (sniax)^ /A^ for any node to 
be coirnected to a spanniirg cluster, even at p = Pc = ^, 
as shown in Fig. |31 In fact, the continuity of P^o (p) is 
interrupted merely because T* = 1 becomes an unsta- 
ble fixed point of Eq. ([IJ below pc- There the stable 
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Figure 4: Depiction of the non-planar Hanoi network HNNP 
of finite size^. Starting from a Id-lattice (horizontal lines), 
a set of long-range bonds (arced lines) is added hierarchi- 
cally that break planarity. RG on this network remains exact. 
Nodes on the lowest two levels of the hierarchy have degree 3, 
then degree 5, 7, etc., comprising a fraction of 1/2, 1/4, 1/8, 
etc., of all nodes, resulting in an average degree 4. 

branch transitions to Eq. ^ that lacks extensive clus- 
ters. Hence, it is the intersection of two stable branches 
of fixed points T* that causes a discontinuous transition. 

We believe, that such a transition is generic for this 
type of hierarchical network. For instance, on the non- 
planar Hanoi network HNNP- described in Fig. IH an 
equivalent but far more elaborate calculation results in a 



pi = 0.319445181 ... as the root of an octic polynomial in 
p, while p„ = (3 - /2 = 0.381966 .... We find = 
below pi , at which the exponent ^ jumps to a finite value, 
followed by a domain with In 1.66639/ ln2 « 0.736726 < 
^ ip) < 1 for pi < p < Pu- At ptj, again, two lines of 
fixed points intersect, resulting in an extensive cluster, 
^I' = 1, above and at pu = Pc, causing a discontinuity in 
Poo- For p ^ Pc, finite-size corrections decay as A^*(p)~i, 
where 1 — (p) ^ [pc — p) with a k > 1 , signaling sub- 
tle variations in the cluster formation dynamics between 
networks. 

We have shown that a hierarchy of small- world bonds 
grafted onto a finitely ramified lattice results generically 
in an explosive percolation transition. This differs signifi- 
cantly from other models, where a non-local set of bonds 
or a global growth parameter is considered to achieve 
explosive percolation^ii^"— . At this point, the precise 
conditions to be imposed on the hierarchy of long-range 
bonds for obtaining this transition are not entirely clear. 
However, they have to result in an infinitely ramified 
network^'^ to provide pc < 1, as the example of the Hanoi 
network HN3 demonstrates^. As well, a finitely ramified 
lattice seems essential, as studies of hierarchical networks 
with small-world bonds imposed on infinitely ramified 2d- 
lattice a -'-ii'""' appear to result in infinite-order transitions, 
which have been observed in many other networks^. 
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On a lattice, a fractal only exists in exactly one point, 

pi =Pu = Pc- 



